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1 .  INTRODUCTION 

The  injection  and  propagation  of  a  charge-neutral  beam  across  a  magnetic 

1~  9 

field  has  been  the  subject  of  several  theoretical  and  experimental  studies. 

Such  beams  have  been  considered  for  use  in  tokamak  heating1'  and  driving  cur¬ 
rents  in  a  magnetically  confined  plasma.11  Charge-neutral  beams  have  also  begun 

12—17 

to  supplement  electron  beams  as  a  means  to  actively  probe  the  ionosphere. 

Our  interest  in  the  problem  of  injecting  a  plasma  beam  across  a  magnetic 
field  stems  from  the  need  to  control  spacecraft  charging  during  electron  beam 
experiments.  A  spacecraft  may  not  be  able  to  collect  sufficient  numbers  of 
electrons  from  the  ionosphere  to  compensate  for  those  emitted  in  the  beam.  A 
plasma  beam  injected  transverse  to  the  ambient  magnetic  field  might  serve  to 
collect  electrons  traveling  along  the  field  lines  during  electron  beam 
experiments.  The  collecting  area  of  the  plasma  beam  would  depend  on  how  far  it 
could  penetrate  the  field.  This  idea  is  closely  allied  to  that  of  the  plasma 
contactor  currently  in  use.  There  is  also  the  possibiity  of  using  a  charge- 
neutral  beam  as  an  ULF  antenna  provided  the  beam  can  be  made  long  enough. 

A  beam  injected  into  a  transverse  magnetic  field  may  be  categorized  by  the 
ratio  of  its  radius  R  to  that  of  an  ion  gyroradius  P^.4’5  The  case  of  R/p^  >  1 
is  termed  a  sraal  1-gyroradius  beam  while  that  of  R/p-[  <  1  is  termed  a  large- 
gyroradius  beam.  The  mechanism  by  which  the  beam  penetrates  the  magnetic  field 
also  serves  to  characterize  it.  The  completely  diamagnetic  plasma  with  a 
kinetic  3  »  1  represents  one  extreme.  For  this  case  the  plasma  expels  the 
magnetic  field  from  its  interior  and  Dropagates  ball ist ical lv . 18  20  The  other 
extreme  is  the  beam  which  is  completely  nond i amagnet ic  as  it  penetrates  the 
transverse  field.  In  the  nond iamagnet ic  case  the  beam  may  propagate  by  means 
of  electric  polarization.18,21’22  The  electric  polarization  of  the  beam  arises 
from  space-charge  separation  layers  forming  due  to  the  Lorentz  force.  The  work 
to  be  presented  here  is  restricted  to  the  nond iaraagne tic  regime  for  the  large- 
gyroradius  case. 

The  beam  geometry  used  is  shown  in  Figure  l.  The  beam  is  injected  in  the 
x  direction  with  velocity  v{,  transverse  to  a  uniform  magnetic  field  Bz.  The 
beam  has  a  finite  height  h  in  the  v  direction  and  is  infinite  along  z.  When  the 
beam  crosses  the  magnetic  field  the  electrons  and  ions  in  the  head  of  the  beam 
initially  gyrate  in  opposite  directions.  The  ensuing  current  Jv  creates  space- 
charge  boundary  layers  on  either  side  of  the  beam.  The  accumulation  of  charge 
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on  Che  boundaries  gives  rise  to  a  transverse  electric  field  Ev .  This  polariza- 
tion  field  will  result  in  the  beam  particles  acquiring  an  E  *  3  drift  across  the 
magnetic  field. 

A  theory  for  the  steady-state  polarization  drift  has  been  determined  by 
22 

Schmidt.  The  first  requirement  is  that  the  beam  possess  a  kinetic  energy 
density  much  greater  than  the  polarization  field  energy  density. 

The  additional  requirement  that  the  beam  E  *  B  drift  with  aoproximately  the 
initial  injection  velocity  yields 

£  -  1  >  »  1  ,  (  1  ) 


where  £  is  the  dielectric  constant  of  the  beam  plasma,  is  the  ion  plasma 

frequency,  and  is  the  ion  gyrof requencv .  It  has  been  further  shown3  that  a 
condition  for  charge  neutrality  is 


l/2 

mb  i 

£  »  (— ) 


mc 
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where  m^^/mp  is  Che  beam  ion  to  electron  mass  ratio. 

A  second  necessary  condition  for  the  polarization  drift  model  to  work  is 
that  the  space-charge  boundary  layer  be  significantly  smaller  than  the  beam 
height.  This  is  expressed  as 


<5  =  p  £_/(  e  —  1)  «  h.  (3) 

This  represents  the  lower  limit  of  the  beam  height.  An  upper  limit  is 

established  by  requiring  that  the  maximum  potential  difference  across  the  beam 

.  .  .  .  23 

cannot  exceed  the  initial  kinetic  energy  which  the  10ns  are  accelerated  to. 

2 

Thus,  eE <  Mv^  / 2 ,  and  may  be  written  as 

h  <  Pi/2.  (4) 
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A  charge-neutral  beam  is  expected  to  propagate  across  a  transverse  field  if 
these  conditions  are  satisfied.  The  time-averaged  drift  velocity  of  the  beam^ 
is  given  by 

vx  =  vb  (1-1/e).  (5) 

The  drift  velocity  is  very  nearly  the  injection  velocity  for  e  »  1. 

In  section  2  the  parameters  used  in  the  simulations  are  described.  Beam 
injection  into  a  vacuum  is  presented  in  section  3.  For  vacuum  injection  two 
cases  are  considered.  The  parameters  of  the  first  case  do  not  satisfy  the 
requirements  for  the  polarization  drift  model  and  the  beam  separates  into  its 
charged  components.  For  che  second  case  the  polarization  requirements  are 
satisfied  and  the  beam  does  propagate.  However,  the  beam  curves  slightly  in  the 
direction  of  initial  electron  gyration.  Section  4  deals  with  injection  into  a 
tenuous  magnetized  plasma  where  the  mass  of  the  plasma  ions  is  equal  to  the  mass 
of  the  beam  ions.  With  an  ambient  plasma  present  of  density  nD/nb  =  1/100,  the 
curvature  of  the  beam  is  noticeably  increased.  For  the  case  of  np/nb  =  1/10, 
the  beam  curves  sharply  in  the  direction  of  electron  gyration  and  then  in  the 
opposite  direction.  Results  are  summarized  in  section  5. 


2.  SIMULATION  PARAMETERS 

The  injection  of  the  beam  and  the  response  of  the  bean  and  ambient  plasma 
are  investigated  with  two-dimensional  (three-velocitv)  electrostatic  full- 
particle  simulations.  The  simulation  model  treats  the  beam  source,  the  beam 
itself,  and  any  ambient  plasma  as  an  isolated  system;  periodic  boundary  condi¬ 
tions  are  not  employed.  The  simulation  code  is  the  same  as  used  by  Pritchett 

.  2  4  2  5 

and  Winglee  for  electron  beam  studies  and  Winglee  and  Pritchett  for  charge- 

25 

neutral  beam  studies.  In  this  study,  as  in  the  study  by  Winglee  and  Pritchett  , 
the  beam  is  charge-neutral  rather  than  consisting  solely  of  electrons.  The  mass 
of  the  beam  ions  is  100  times  the  electron  mass  me.  When  ambient  plasma  is 

present  the  mass  of  the  beam  ions  m^i  is  equal  to  the  mass  of  the  ambient  ions 
tripi  .  The  magnetic  field  is  directed  along  the  z  axis  (out  of  the  simulation 
plane)  and  is  constant  throughout  the  volume  of  the  simulation  box.  In  most  of 
the  examples  presented  the  magnetic  field  strength  is  such  that  the  ratio  of 
electron  cyclotron  frequency  12e  to  the  beam-electron  plasma  frequency  m^g  is 

=  1/4.  The  unit  of  length  in  the  simulations  is  the  grid  spacing  A  and 
the  unit  of  time  is  the  inverse  beam-electron  plasma  frequency  where 

•°be  =  (nbOe  /e0mbe)  >  and  nb0  *-s  c^e  beam  density  at  injection.  The  initial 
velocity  v],  is  typically  0.5A<Jbe.  An  ion  gyroradius  is  typically  200A.  For  a 
400  eV  Argon  beam  injected  across  a  0.25  gauss  field  A  is  about  3.6  m.  The 
system  dimensions  are  Lx  =  512A  and  Ly  =  64A  or  128A.  The  time  step  in  the 
simulation  is  typically  At  =  0.1/ojjjg.  The  largest  system  of  particles  consists 
of  224,000  particles  and  the  smallest  consists  of  138,500.  The  beam  at  injec¬ 
tion  has  a  height  h  =  10A  perpendicular  to  the  magnetic  field.  For  2e  /  ui^e  =  1/4 
the  beam  height  is  20  times  the  beam  electron  gyroradius  and  1/20  of  the  beam 
ion  gyroradius  for  the  parameters  considered  in  this  work.  The  thermal  electron 
velocity  is  typically  0.5A(jj3e  and  the  thermal  ion  velocity  is  typically 

0 . 05AL013P  . 

The  beam  is  injected  from  the  center  of  the  right-hand  side  of  an  idealized 
spacecraft.  A  particle  which  encounters  the  boundary  of  the  spacecraft  is 
removed  from  the  svscem  and  its  charge  uniformly  distributed  on  the  surface  of 
the  spacecraft.  The  beam  is  taken  to  be  charge-neutral  so  the  injection  process 
leaves  no  net  charge  on  the  spacecraft.  However,  spacecraft  charging  can  occur 
if  particles  are  accelerated  onto  the  spacecraft  by  fields  produced  by 
space-charge  effects.  The  spacecraft  has  the  dimensions  of  6A  *  12A  to  allow 
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for  the  collection  of  such  a  charge.  The  center  of  the  spacecraft  is  located  at 
(0.11  Lx,  0.50  Ly). 

Ambient  plasma  is  reflected  from  the  boundaries  of  the  simulation  box  while 
beam  particles  are  allowed  to  escape.  This  feature  of  the  simulation  model 
remains  valid  provided  the  plasma  is  basically  undisturbed  near  the  boundary.  A 
loss  of  4 7*  of  the  beam  particles  is  considered  excessive  and  the  simulations  are 
assumed  invalid  beyond  that  point. 
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3.  INJECTION  INTO  VACUUM 

In  this  section  we  examine  the  behavior  of  a  charge-neutral  beam  injected 
transverse  to  a  constant  vacuum  magnetic  field.  Two  cases  are  considered.  The 
first  has  parameters  which  violate  some  of  the  polarization  drift  model 
criteria.  As  ^uch,  it  would  not  be  expected  to  propagate  through  the  magnetic 
field  by  virtue  of  collective  effects.  The  second  case  is  typical  of  the 
simulation  parameters  presented  in  section  2  and  does  exhibit  the  expected 
polarization  drift. 

3.1  Nonpropagating  3eam.  For  the  nonpropagating  example  the  injection 

velocity  v^,  =  lOAoa^g  and  the  magnetic  field  is  =  1.0.  The  beam  height 

is  h  =  10A  .  The  beam  ion  mass  to  electron  mass  ratio  is  mv,i  /me  =  100.  The  ion 
gyroradius  is  then  =  1000A  .  The  dielectric  constant  of  the  beam  is  e  =  100. 
This  marginally  satisfies  the  requirement  e  »  The  space-charge 

boundary  layer  thickness  is  5  =  10A  and  is  equal  to  the  beam  height,  a  clear 
violation  of  the  polarization  model  criteria. 

Figure  2  shows  contour  plots  of  charge  density  for  this  case  shortly  after 
injection  is  initiated.  The  dotted  contours  indicate  regions  where  electron 
charge  density  dominates  while  the  solid  contours  indicate  regions  dominated  by 
positive  charge.  The  spacecraft  is  negatively  charged  and  represented  by  the 
rectangular  contours  behind  the  injection  plane  at  x  =  60A.  Upon  injection  the 
initially  charge-neutral  beam  separates  into  charged  beams,  an  ion  beam  and  an 
electron  beam,  which  gyrate  in  the  magnetic  field.  The  beam  polarization  drift 
is  not  observed  in  this  case. 

3.2  Propagating  Beam.  The  example  of  the  propagating  beam  is  typical  of 
the  simulation  parameters  used  throughout  the  rest  of  this  study.  The 
dielectric  constant  for  this  case  is  e  =  1600  which  easily  satisfies  the 

requirement  £  >  (mv^/mg)1  =  10.  The  space-charge  boundary  layer  thickness  is 

5  =  0.125A.  The  beam  height  is  h  =  10A ;  consequently  the  requirement  that  5  «  h 
is  satisfied. 

3.2.1  Beam  density  and  potential  contour  plots.  Figure  3  shows  separate 
beam  density  contour  plots  for  the  electrons  and  ions.  The  minimum  contour 
value  in  these  plots  is  10%  of  the  injection  density.  During  a  period  of 
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Electron  Density  Contours 
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Figure  3 


Example  of  beam  propagating  in  vacuum.  Minimum  density 
contour  is  10%  of  injection  value.  Beam  curves  upward 
s 1 ightlv . 


880wjje”  the  head  of  the  beam  travels  from  the  initial  injection  plane  at 
x  *  6QA  to  approximately  x  **  500A .  This  yields  an  average  propagation  velocity 
of  vx  *  0.5Au>tje,  which  is  the  initial  injection  velocity,  an  expected  result. 
Close  examination  of  the  beam  contours  reveals  a  slight  curvature  of  the  beam 
path  of  both  ions  and  electrons.  The  beam  is  deflected  in  the  y  direction, 
which  is  opposite  to  that  in  which  the  beam  ions  would  initially  gyrate. 

At  the  end  of  the  run  the  head  of  the  beam  has  expanded  to  a  height  of 
approximately  h  =  15A  and  traveled  x  =  2.2p{.  Both  spreading  and  particle  loss 
from  the  space-charge  layers  degrade  the  density  of  the  charge-netural  region 
within  the  beam.  The  density  of  the  charge-neutral  region  has  been  estimated  at 
x  =  161A,  261A,  and  361A  and  shows  a  noticeable  decrease.  Linear  extrapolation 
of  these  density  values  suggests  the  beam  will  degrade  to  10%  of  the  injection 
value  after  traveling  4.5p^. 

Figure  4  shows  contour  plots  of  the  charge  density  with  the  minimum  contour 
line  equal  to  1%  of  the  injection  density.  As  before,  the  dotted  lines  indicate 
regions  of  negative  charge  density  and  the  solid  lines  indicate  positive  charge 
density  regions.  The  beam  is  shown  evolving  over  a  period  of  160cu^e-1.  The 
most  striking  feature  of  this  plot  is  the  asymmetric  loss  of  ions  from  the 
lower,  or  positively  charged,  side  of  the  beam.  Ions  at  the  outer  surface  of 
the  positive  space-charge  boundary  layer  are  partially  shielded  from  the 
electric  polarization  field  and  can  eventually  gyrate  away  from  the  beam.  The 
dashed  line  on  the  lower  side  of  the  beam  represents  a  gyro-orbit  of  an  ion 
moving  at  the  injection  velocity.  An  additional  noteworthy  feature  is  the 
conspicuous  bulge  of  electrons  above  the  head  of  the  beam. 

Electric  potential  contour  plots  to  accompany  Figure  4  are  shown  in  Figure 
o.  As  the  beam  transverses  the  magnetic  field  a  potential  distribution  forms 
which  has  a  gradient  perpend icul ar  to  both  the  magnetic  field  and  the  direction 
of  beam  injection.  The  dotted  lines  indicate  regions  of  negative  potential  and 
the  solid  lines  indicate  regions  of  positive  potential.  From  the  first  two 
panels  of  the  plot  it  is  evident  that  within  the  beam  region  the  electric  field 
is  primarily  in  the  y  direction.  The  gradient  in  the  potential  contours  at  the 
front,  or  head,  of  the  beam  indicate  a  rearward  pointing  electric  field  above 
the  beam  and  a  forward  pointing  field  below  the  beam.  The  electrons  on  the 
outer  edge  of  the  beam  head  are  extruded  outward  from  the  beam  by  virtue  of  the 
resulting  E  *  3  drift  and  produces  the  electron  bulge. 
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sequence  of  potential  lines  for  the  vacuum 


The  asymmetric  ion  loss  from  the  beam  serves  to  distort  the  potential 
distribution.  As  these  ions  gyrate  away  from  the  beam  the  polarization  layer  is 
broadened.  The  polarization  electric  field  then  extends  over  a  region  much 
higher  than  the  original  beam  height  at  injection. 

3.2.2  3eam  density  and  current  plots.  In  addition  to  distorting  the 
potential  distribution,  the  loss  of  ions  from  the  space-charge  layer  allows  a 
net  current  to  flow  within  the  beam.  The  beam  electrons  drifting  across  the 
magnetic  field  constitute  a  current  flowing  from  the  head  of  the  beam  toward  the 

4  *► 

source.  The  sense  of  the  resulting  j  *  B  force  is  in  the  correct  direction  to 
account  for  the  transverse  displacement.  This  is  demonstrated  in  Figure  6. 

Here  the  10%  density  contour  of  the  ion  component  is  shown  below  the  net 
current  flow  in  the  beam.  The  forward  half  of  the  beam  is  discernibly  shifted 
in  the  y  direction.  Within  this  portion  of  the  beam  a  net  current  in  the  -x 
direction  is  evident.  Additionally,  a  current  structure  flowing  from  the  head 
to  the  source  persists  along  the  entire  length  of  the  positively  charged  side  of 
the  beam. 
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4.  INJECTION  INTO  A  LOW-DENSITY  PLASMA 

The  presence  of  an  ambient  plasma  provides  a  source  of  particles  which 
responds  to  the  locally  induced  fields  of  the  beam  plasma.  The  resulting  flow 
of  the  ambient  particles  serves  to  maintain  quasi-neutrality  in  the  system.  The 
polarization  drift  model  relies  on  charge  imbalance  to  drive  the  beam  across  the 
magnetic  field.  The  ambient  plasma  can  partially  neutralize  the  space-charge 
boundary  layers,  altering  both  the  trajectory  and  profile  of  the  beam. 

Injection  as  a  function  of  the  ratio  of  ambient  to  beam  plasma  density  is  con¬ 
sidered  in  this  section.  The  cases  of  np/n^  =  1/100  and  np/njj  =  1/10  are  pre¬ 
sented  and  contrasted  with  the  vacuum  injection  case. 

4.1  Beam  Properties  as  a  Function  of  Relative  Density 

4.1.1  Beam  density  and  potential  contour  plots  for  ambient  plasma.  Figure 
7  depicts  the  beam  contour  plots  for  the  case  of.np/n^,  =  1/100.  As  in  Figure  3 
the  minimum  contour  value  is  10%  of  the  injection  value.  During  a  period  of 
720wbe_1  the  beam  travels  from  x  3  604  to  x  =  4204.  As  in  the  vacuum  case,  the 
average  propagation  velocity  is  vx  =  0.5wbe,  the  injection  velocity.  In  this 
case  the  beam  curvature  is  visibly  more  pronounced  than  the  vacuum  case. 

Density  estimates  along  the  beam  suggest  it  will  degrade  to  10%  of  tne  injection 
value  after  traveling  2.7p^. 

The  density  contour  plots  for  the  Op/n^  =  1/10  case  are  presented  in 
Figure  3.  Here  there  is  a  gross  distortion  of  the  expected  beam  trajectory. 
Close  inspection  of  the  two  panels  shows  the  beam  species  to  literally  separate 
as  the  beam  initially  deflects  in  the  v  direction.  The  two  streams  then 
recombine  and  the  beam  deflects  in  the  opposite  direction.  However,  the  head  of 

a  1 

the  beam  still  has  an  average  velocity  in  the  x  direction  of  vx  =  0.54ojbe-1.  No 
estimate  is  made  of  the  distance  this  beam  might  travel. 

Beam  density  degradation  as  a  function  of  ambient  plasma  density  is  presen¬ 
ted  in  Figure  9.  The  density  profiles  of  the  beam  cross-section  are  presented 
in  order  of  increasing  ambient  density  with  the  vacuum  case  at  the  ton.  In  all 
cases  the  beam  has  propagated  for  a  time  of  320wv)e~  1  and  the  cross-sect  ion  is 
taken  at  x  =  1614.  The  beam  is  propagating  out  of  the  page  in  Figure  9.  The 
vertical  axis  is  scaled  to  the  normalized  beam  densitv  and  the  horizontal  axis 
is  the  y  axis  of  the  simulation  box.  The  dashed  line  forming  a  rectangle 
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Figure  7.  Example  of  beam  propagating  in  tenuous  plasma  (rtp/n^,  = 
1/100).  Minimum  density  contour  is  10%  of  injection 
value.  Beam  curvature  is  more  pronounced  than  Che  vacuum 
case  . 
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Comparison  of  density  profiles  for  beams  which  have 

propagated  for  t  3  320  Cross-sections  are  taken  at 

x  J  161  A.  (a)  Vacuum,  (b)  Tenuous  plasma,  (c)  Marginally 
dense  plasma. 
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represents  the  idealized  density  profile  of  the  beam  cross-section  at  injec¬ 
tion.  The  solid  line  represents  the  density  profile  of  the  beam  ions.  The 
dotted  line  represents  the  density  profile  of  the  beam  electrons.  For  the 
vacuum  case  shown  in  panel  (a)  the  beam  is  slightly  offset  to  the  right  of  the 
initial  injection  profile.  Ions  in  noticeable  quantities  are  moving  to  the  left 
of  the  positive  side  of  the  beam  to  a  distance  of  approximately  20A.  For  the 
tenuous  plasma  case  shown  in  panel  (b)  the  beam  is  offset  to  a  greater  degree 
than  the  vacuum  case.  Also,  ions  moving  to  the  left  of  the  positive  side  of  the 
beam  extend  beyond  the  boundary  of  the  simulation  box.  In  the  case  of  the 
marginally  dense  plasma  shown  in  panel  (c)  the  beam  is  almost  completely  outside 
of  the  initial  injection  profile.  The  ions  moving  to  the  left  do  so  in  greater 
quantity  than  the  previous  two  cases  and  extend  Keyond  the  boundary  of  the 
simulation  box.  In  all  three  cases  fluctuations  in  the  ion  density  are  followed 
closely  by  electron  density  fluctuations. 

The  increase  of  beam  curvature  for  the  cases  of  interpenetrating  ambient 
plasma  is  clarified  in  Figure  10.  In  this  example  the  10%  density  contour  is 
shown  below  the  net  beam  current  for  Op/n^,  =  1/100.  The  beam  has  been  propa¬ 
gating  for  t  *  320  o)be_1<  A  backward  flowing  current  structure  has  developed 
along  the  bottom  of  the  beam.  The  current  is  roughly  twice  the  magnitude  of  the 
vacuum  case  shown  in  Figure  6  and  has  evolved  in  half  the  time.  The  j  *  B 
force  is  therefore  larger  and  can  act  on  the  beam  for  a  longer  time. 

4.1.2  Potential  contour  plots.  Potential  distortion  as  a  function. of 
ambient  plasma  density  is  presented  in  Figure  11.  The  potential  plots  depicted 
are  companions  to  those  in  Figure  9.  The  vacuum  case  occupies  the  left  panel. 
The  np/nb  =  1/100  case  in  the  center  panel  shows  the  potential  contours  have 
been  distended  downward  in  the  direction  in  which  ions  are  being  shed  from  the 
positive  space-charge  boundary  layer.  The  contours  are  tilted  slightly  upward, 
coinciding  with  the  beam  trajectory.  The  right  panel  depicts  the  np/n),  =  1/10 
case.  The  potential  contours  are  grossly  distorted  when  compared  to  the  vacuum 
case.  The  beam  trajectory,  however,  is  still  along  the  contours  which  indicates 
an  E  x  3  drift  mechanism. 

4.2  Response  of  the  Ambient  Plasma.  From  study  of  the  vacuum  case  it  is 
obvious  that  the  shedding  of  ions  is  partially  responsible  for  the  potential 
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Figure  10.  Beam  current  plot  for  tenuous  plasma  case  (n^/nh  =  1/100). 

Ion  density  contour  and  current  plots  are  shown  separately. 

•f 

Current  in  lower  edge  appears  Co  produce  j  *  B  force  to 
curve  beam  upward. 
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beams  which  nave  propagated  for 


contour  distortion.  The  presence  of  ambient  plasma  enhances  this  beam  ion  loss 
by  providing  additional  shielding  of  the  space-charge  layers.  This  is  evident 
from  a  consideration  of  Figure  12.  The  figure  depicts  the  ambient  currents  for 
the  np/n^  -  1/100  case.  The  current  of  each  species  is  shown  separately.  The 
beam  region  is  identified  by  the  superimposed  density  contour.  Ambient  ion 
currents  are  shown  in  the  upper  panel  and  ambient  electron  currents  in  the  lower 
panel.  The  obvious  difference  of  the  two  plots  is  attributable  to  magnetization 
scale  lengths.  For  distances  on  the  order  of  the  beam  height  the  ions  are 
virtually  unmagnetized  whereas  the  electrons  are  strongly  magnetized.  When 
acted  on  by  an  electric  field  the  ions  undergo  direct  acceleration  while  the 

electrons  1  x  tj  drift. 

Within  the  beam  region  ambient  ions  are  driven  by  the  polarization  field. 
This  accounts  for  the  ion  current  observed  flowing  across  the  beam  in  Figure 
12.  The  ambient  ion  flow  results  in  the  transverse  polarization  field  being 
reduced.  Below  the  beam  the  ambient  ion  current  "is  driven  by  the  field  of  the 
beam  ions  being  shed  from  the  positive  space-charge  layer.  The  ambient  ions  are 
also  drawn  to  the  spacecraft  since  it  has  become  negatively  charged. 

On  the  lower  side  of  the  beam  the  ambient  electron  current  ’is  evidence  of 
the  polarization  region  being  distended  downward.  The  electrons  drifting  in  the 
x  direction  imply  an  electric  field  in  the  y  direction. 

The  absence  of  ambient  currents  in  much  of  the  beam  region  is  due  to  the 
expulsion  of  the  ambient  plasma  by  the  mechanisms  described  above. 
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Figure  12.  Ambient  current  plot  for  tenuous  plasma  case  (nn/nh  =1/100). 


Ion  density  contour  is  superimposed,  (a)  Ambient  ion  current 
(b)  Ambient  electron  current. 


5.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  investigated  the  properties  of  a  charge-neutral  beam 
injected  into  a  transverse  magnetic  field.  The  beam  consisted  of  ions  and 
electrons.  The  study  made  use  of  two-dimensional  electrostatic  simulations 
employing  an  isolated  system  model.  Motivation  for  the  study  stems  from  an 
interest  in  creating  long  plasma  columns  in  the  ionosphere.  For  a  limited  range 
of  parameters  it  is  known  that  a  charge-neutral  beam  can  propagate  perpendicular 
to  a  magnetic  field  by  virtue  of  collective  effects.  We  have  simulated  this 
phenomenon  to  determine  possible  limitations  on  the  distance  a  beam  might  travel 
across  a  magnetic  field.  It  has  been  shown  that  when  the  criteria  for  the 
polarization  drift  model  are  satisfied,  the  beam  does  propagate. 

For  the  case  of  vacuum  injection  (op/n^  =  0)  the  beam  propagates  at  the 
injection  velocity.  The  density  of  the  charge-neutral  core  of  the  beam  is 
depleted  by  a  combination  of  beam  spreading  and  particle  loss.  It  is  estimated 
that  the  beam  core  will  degrade  to  10%  of  its  injection  density  after  traveling 
4.5p^.  For  a  400  eV  Argon  beam  injected  across  a  0.25  gauss  field  this  is  more 
than  3  kilometers.  The  beam  exhibits  a  slight  curvature  in  the  direction 
opposite  to  that  of  ion  gyration  in  the  magnetic  field.  The  large  gyroradius  of 
the  ions  allow  some  of  the  outer  edge  of  the  beam  to  be  lost  from  the  positive 
space-charge  layer.  This  leaves  a  net  electron  current  in  the  beam  flowing 
toward  the  source  from  the  head.  The  resulting  j  *  B  force  is  in  the  correct 
direction  to  account  for  the  beam  curvature. 

In  the  case  of  injection  into  a  tenuous  plasma  (np/n^  =  1/100)  the  beam 

still  propagates  at  the  injection  velocity.  The  beam  is  estimated  to  degrade  to 

10%  of  the  injection  density  after  traveling  2.7p^.  Beam  curvature  is  again 

exhibited  and  is  more  pronounced  than  for  the  vacuum  case.  The  electron  current 

in  the  beam  is  greater  than  in  the  vacuum  case  and  is  established  sooner.  The 
+  *► 

j  x  B  force  would  then  be  larger,  accounting  for  the  increased  curving  of  the 
beam.  The  increase  in  beam  electron  current  is  due  to  the  enhanced  shielding 
provided  h»y  the  ambient  plasma.  The  ambient  ions  flow  directly  across  the  beam, 
partially  shorting  the  polarization  field.  This  aids  the  loss  of  positive 
charge  from  the  lower  space-charge  polarization  layer.  The  escaping  ions 
noticeably  distend  the.  potential  contours  and  draw  the  polarization  field  down 
from  the  beam. 
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Injection  into  a  marginally  dense  (np/nb  =  1/10)  plasma  yields  a  propaga¬ 
tion  speed  in  the  x  direction  equal  to  the  injection  velocity.  No  estimate  of 
propagation  distance  was  made  for  this  case.  The  beam  curves  drastically  and 
partially  separates  into  distinct  streams  of  ions  and  electrons.  The  loss  of 
ions  from  the  beam  exceeds  that  of  the  tenuous  case  and  the  potential  contours 
deviate  markedly  from  those  of  the  other  cases. 

The  presence  of  the  expected  polarization  field  required  to  drive  the  beam 
through  a  magnetic  barrier  has  been  confirmed  by  our  simulations.  The  beam 
degradation  reported  by  other  investigators^  has  also  been  confirmed.  The  net 
electron  current  in  the  beam  and  the  subsequent  beam  curvature  were  unexpected. 

Despite  the  rapid  degradation  of  the  beam  it  appears  possible  to  produce 
extremely  long  plasma  columns  using  the  transverse  injection  technique  studied 
in  this  work.  A  test  of  the  feasibility  of  using  such  a  column  as  an  electron 
collector  would  require  a  full  three-dimensional  simulation. 
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MISSION 

°f 

Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (C*I)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  CSI  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


